
% Annual
m=nanmean(nmodel.data(:,:,1:12),3);
d=nanmean(dat.data(:,:,1:12),3);
di=d-m;
%load cmap_insitu
close all
gccolor('mp',dat.londata,dat.latdata,m','on',0,1000,['Model - Annual - ',year,' (mgC.m-2.d-1) '])
printgc('-f1','-depsc2',['./PLOT_VALIDATION/PP/NPP_model_',parameter,'_',bbio,'_y',year,'_ANNUAL.eps'])
%colormap(cmap_insitu)
close all
gccolor('mp',dat.londata,dat.latdata,d','on',0,1000,['Data - Annual - ',year,' (mgC.m-2.d-1) '])
printgc('-f1','-depsc2',['./PLOT_VALIDATION/PP/NPP_data_',parameter,'_',bbio,'_y',year,'_ANNUAL.eps'])
%colormap(cmap_insitu)
load cmap_politoliu
close all
gccolor('mp',dat.londata,dat.latdata,di','on',-1000,1000,['(data-model) - Annual - ',year,' (mgC.m-2.d-1) '])
colormap(politoliu)
printgc('-f1','-depsc2',['./PLOT_VALIDATION/PP/NPP_diff_',parameter,'_',bbio,'_y',year,'_ANNUAL.eps'])

% JFM
m=nanmean(nmodel.data(:,:,1:3),3);
d=nanmean(dat.data(:,:,1:3),3);
di=d-m;
%load cmap_insitu
close all
gccolor('mp',dat.londata,dat.latdata,m','on',0,1000,['Model - JFM - ',year,' (mgC.m-2.d-1) '])
printgc('-f1','-depsc2',['./PLOT_VALIDATION/PP/NPP_model_',parameter,'_',bbio,'_y',year,'_JFM.eps'])
%colormap(cmap_insitu)
close all
gccolor('mp',dat.londata,dat.latdata,d','on',0,1000,['Data - JFM - ',year,' (mgC.m-2.d-1) '])
printgc('-f1','-depsc2',['./PLOT_VALIDATION/PP/NPP_data_',parameter,'_',bbio,'_y',year,'_JFM.eps'])
%colormap(cmap_insitu)
load cmap_politoliu
close all
gccolor('mp',dat.londata,dat.latdata,di','on',-1000,1000,['(data-model) - JFM - ',year,' (mgC.m-2.d-1) '])
colormap(politoliu)
printgc('-f1','-depsc2',['./PLOT_VALIDATION/PP/NPP_diff_',parameter,'_',bbio,'_y',year,'_JFM.eps'])

% AMJ
m=nanmean(nmodel.data(:,:,4:6),3);
d=nanmean(dat.data(:,:,4:6),3);
di=d-m;
%load cmap_insitu
close all
gccolor('mp',dat.londata,dat.latdata,m','on',0,1000,['Model - AMJ - ',year,' (mgC.m-2.d-1) '])
printgc('-f1','-depsc2',['./PLOT_VALIDATION/PP/NPP_model_',parameter,'_',bbio,'_y',year,'_AMJ.eps'])
%colormap(cmap_insitu)
close all
gccolor('mp',dat.londata,dat.latdata,d','on',0,1000,['Data - AMJ - ',year,' (mgC.m-2.d-1) '])
printgc('-f1','-depsc2',['./PLOT_VALIDATION/PP/NPP_data_',parameter,'_',bbio,'_y',year,'_AMJ.eps'])
%colormap(cmap_insitu)
load cmap_politoliu
close all
gccolor('mp',dat.londata,dat.latdata,di','on',-1000,1000,['(data-model) - AMJ - ',year,' (mgC.m-2.d-1) '])
colormap(politoliu)
printgc('-f1','-depsc2',['./PLOT_VALIDATION/PP/NPP_diff_',parameter,'_',bbio,'_y',year,'_AMJ.eps'])

% JAS
m=nanmean(nmodel.data(:,:,7:9),3);
d=nanmean(dat.data(:,:,7:9),3);
di=d-m;
%load cmap_insitu
close all
gccolor('mp',dat.londata,dat.latdata,m','on',0,1000,['Model - JAS - ',year,' (mgC.m-2.d-1) '])
printgc('-f1','-depsc2',['./PLOT_VALIDATION/PP/NPP_model_',parameter,'_',bbio,'_y',year,'_JAS.eps'])
%colormap(cmap_insitu)
close all
gccolor('mp',dat.londata,dat.latdata,d','on',0,1000,['Data - JAS - ',year,' (mgC.m-2.d-1) '])
printgc('-f1','-depsc2',['./PLOT_VALIDATION/PP/NPP_data_',parameter,'_',bbio,'_y',year,'_JAS.eps'])
%colormap(cmap_insitu)
load cmap_politoliu
close all
gccolor('mp',dat.londata,dat.latdata,di','on',-1000,1000,['(data-model) - JAS - ',year,' (mgC.m-2.d-1) '])
colormap(politoliu)
printgc('-f1','-depsc2',['./PLOT_VALIDATION/PP/NPP_diff_',parameter,'_',bbio,'_y',year,'_JAS.eps'])

% OND
m=nanmean(nmodel.data(:,:,10:12),3);
d=nanmean(dat.data(:,:,10:12),3);
di=d-m;
%load cmap_insitu
close all
gccolor('mp',dat.londata,dat.latdata,m','on',0,1000,['Model - OND - ',year,' (mgC.m-2.d-1) '])
printgc('-f1','-depsc2',['./PLOT_VALIDATION/PP/NPP_model_',parameter,'_',bbio,'_y',year,'_OND.eps'])
%colormap(cmap_insitu)
close all
gccolor('mp',dat.londata,dat.latdata,d','on',0,1000,['Data - OND - ',year,' (mgC.m-2.d-1) '])
printgc('-f1','-depsc2',['./PLOT_VALIDATION/PP/NPP_data_',parameter,'_',bbio,'_y',year,'_OND.eps'])
%colormap(cmap_insitu)
load cmap_politoliu
close all
gccolor('mp',dat.londata,dat.latdata,di','on',-1000,1000,['(data-model) - OND - ',year,' (mgC.m-2.d-1) '])
colormap(politoliu)
printgc('-f1','-depsc2',['./PLOT_VALIDATION/PP/NPP_diff_',parameter,'_',bbio,'_y',year,'_OND.eps'])
close all
